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Abstract. The atmospheric neutrino-induced upward muon 
flux are calculated by using the multi-component primary en- 
ergy spectrum, CORSIKA EAS simulation code for the re- 
production of the atmospheric neutrino spectra and improved 
parton model for charged-current cross sections. The results 
are obtained at 10^ — 10^ GeV muon energy range and 0—89° 
zenith angular range. 



1 Introduction 

Neutrino-induced upward muon flux are both the direct trans- 
mitters of information about neutrino flux from energetic as- 
trophysical sources (binary stars, AGN) and direct posterity 
of atmospheric very high energy neutrino. These problems 
are being investigated in many modern experiments (see re- 
view (Gaisser et al., 1995)) and as a rule, the upward muon 
flux of atmospheric origin is considered in the capacity of the 
background flux. 

Here, on the basis of multi-component model of primary cos- 
mic ray (Biermann, 1993), QGSJET model of high energy 
A — AAir interactions (Kalmykov and Ostapchenko, 1993) 
and parton model predictions of charged-current cross sec- 
tions (Quigg et al., 1986) we evaluated the expected back- 
ground upward muon energy spectra at different zenith an- 
gles and studied (E'^ — E^ n) and (Ef^ — En) correlations. 

2 Method of calculations 

Upward muon differential energy spectra close to Earth sur- 
face at different zenith angles (6) can be represented as fol- 
lows: 



dF{e) 
BE,, 



E 



dE^ 



dEA 



Ea 



5 = 1/^,17^ "'-^f* 



dEf 



■dE, 



dG^ dW 
dE^ dt 



dt 



(1) 




Correspondence to: samvel@jerewanl.yerphi.am 



10" 10 10 
Primary energy (TeV) 

Fig. 1. Primary energy spectra according to 2-component origin of 
cosmic ray. The hatch area contains (Ter-Antonyan and Haroyan, 
2000), DICE (Swordy et al., 2000), and CASA-BLANCA (Fowler 
et al., 2000) data. 



where: 

d'^/EA - is the primary energy spectrum of the nucleus A; 
df{EA, A, 9)/dE^ is the atmospheric neutrino (^ = v^) and 
anti-neutrino (^ = V^) energy spectrum at the primary nu- 
cleus A with Ea energy and zenith angle; 
dG^{t, E^)/ dE^ is the neutrino-induced muon energy spec- 
trum at observation level (close to the Earth surface) at the 
given neutrino energy Ey j: and depth t of muon production 
in the Earth; 

dW{E^^ E^)/dt\s a. probability of charged-current interac- 
tion at t depth of the Earth matter. 

Energy spectra of primary nuclei ( A = \ — 59) according 
to the multi-component model of primary cosmic ray origin 
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Fig. 2. Atmospheric neutrino (symbols) and anti-neutrino (lines) energy spectra at different primary proton energy (Ep,TeV) at zenith 
angles 6» = 0° (left) and 61 = 89° (right). 



(Biermaim, 1993) are presented in 2-component form (see 
Fig.l): 



dEA 



dEA 



■dEA J 



(2) 



where the first component (ISM) is derived from the explo- 
sions of normal supernova into an interstellar medium with 
expected rigidity-dependent power law spectra 
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and the second component (SW) is a result of the explosions 
of stars into their former stellar winds with expected rigidity- 
dependent power law spectra 



^ = / E^^^ ■ Ea < Esw 
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where ^a is a scale factor {E in TeV units) from approxima- 
tions (Wiebel-Sooth and Biermann, 1998), 
EisM = RiSM ■ Z and Egw = ELsw ■ Z are the cor- 
responding rigidity-dependent critical energies (Biermann, 
1993) and Z is the charge of nucleus A. 
The values of spectral parameters from (2-4) are: 



71 = 2.78, 72 = 2.65, 73 = 3.28, Rism ^ 200 TV, 
Rsw — 2000 TV and the fraction of each component: 
Sa,i = 1 - dA,2, 5a,2 = {2ZRisMp^-^^-''^ at 7^=1,0 = 
2.75 and 7A>i,o = 2.66 (Wiebel-Sooth and Biermann, 1998) 
is obtained by normalization of (2-4) with approximation of 
balloon and satellite data (Wiebel-Sooth and Biermann, 1998) 
atEA = l TeV and minimization of KASCADE (Classtet- 
ter et al., 1999) and ANI (ChiHngaryan et al., 1999) EAS size 
spectra with the expected ones using 2-component represen- 
tation (2-4) (Ter-Antonyan and Biermann, 2001). 
Values of primary energy limits in the expression (1) are cho- 
sen Emin = and Emax = 10^ GeV The expected pri- 
mary energy spectra for different nuclei, all-particle and nu- 
cleon energy spectra by 2-component model are shown in 
the Fig. 1. A contribution of the third primary extra-galactic 
component is ignored. 

Atmospheric neutrino (anti-neutrino) energy spectra 
df{EA,A,e)/dE^ at given six A group (A = 1, 4, 12, 
16, 28, 56) and Ea = 3.16 • 10^, lG-\ . . . lO'^ GeV, = 
0, 30, 60, 89° parameters of primary nuclei were calculated 
in the framework of QGSJET interaction model (Kalmykov 
and Ostapchenko, 1993) by CORSIKA562 EAS simulation 
code (Heck et al., 1998). Intermediate values of tabulated 
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spectral data are evaluated by interpolation with preUminary 
linearizations. It should be noted that the neutrino energy 
spectra depend only on primary nucleon energy {Ea/A) with 
accuracy < 5% in A = 1 — 56 range. 

The atmospheric neutrino (anti-neutrino) energy spectra at 
different primary proton energies {Ep = 3.16 • 10^ — 10^ 
GeV) and zenith angles {9 = 0° , 89°) are presented in Fig. 2. 
Neutrino-induced muon energy spectra are determined as 



dEij, a^{y < j/max) dy E^ 
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The probabiUty of charged-current interaction in expression 
(1) at given depth t of the Earth matter has a form: 



dW 

= a^NAexp{-a^NAt) 



(10) 



where cr^ = cTu.yiy < ymax) is a total charged-current cross 
section (7) and Na is the Avogadro number. 
The depth limits of effective (y < ymax) charged-current 
interactions in the Earth matter (t) are 



tmax(6') = / p{r{z))dz. 



(11) 

(12) 



Here p{g/crn?) = 13.8 — 1.73 • 10~^r(cm) is our approxi- 
mation of the Earth density at the radius r, ^max = 2-R cos 6, 
R = 6.371 • 10*cm is the Earth radius, 



and 5^crj/9a;9j/ are corresponding inclusive charged-current L = {l/b)\D.[{a + bE{)/{a + bE^)] 



(13) 



cross sections for the reaction v^.{p^) + iV — > p'^{p~)+ 
anything. Here^ = v^,v^,x = Q"^ /2Mv wdy = v/E^ are 
the scaling variables, is an invariant momentum transfer 
between an incident neutrino and an outgoing muon, v = 
E^ - E* is the energy loss in a laboratory frame, M is the ^ Results 
nucleon mass, Xmin = 5 • 10~^, 



is a rock thickness for muon detection with energy Ef^ on the 
Earth surface at y = 0. 



Umax = 1 - E^/E^, 

E*^ = {a + hEi,) exp[6(t^ax - t)yh - a/h 



(8) 
(9) 



is the muon energy in production point t, parameters a = 

0.002 MeV-cm^/g and 6 = 4 • 10"^ cm^/g (Gaisser et al., 
1995) determine muon energy losses: —dE^/dt = a + bE^. 



Calculations of upward muon energy spectra according to ex- 
pression (1) are performed using numerical integration. The 
charged-current cross section for muon production by neu- 
trino was taken from the renormalization-group-improved par- 
ton model (Quigg et al., 1986) with quark structure functions 
for X > 10""' from (Eichen et al., 1984). We limited the 
range of a; < 10~"^ by Xmin = 5 ■ 10~^ because the check of 
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Fig. 5. (Ei^- < En >) and < E^^j^ >) correlations. 



valence quark structure functions at very small x values by 
quark number sum rules 



/ 
/ 



dx Uv{x, Q^) = 2 



dx dy{x, Q^) = 1 



(14) 



(15) 



showed that the obtained errors of sum rules became more 

than 1 — 2% at x < Xmin- 

Upward muon differential energy spectra at different zenith 
angles are presented in Fig. 3. An accuracy of multidimen- 
sional integration (1) was less than 3%. Corresponding up- 
ward muon integral energy spectra (with accuracy 4 — 5%) 
in comparison with the same calculations (Minorikawa et al., 
1995) are presented in Fig. 4. 



Correlations of muon energy at observation level (near the 
Earth surface) with primary nucleon and parent-neutrino en- 
ergies are shown in Fig. 5. Some irregularities of presented 
data at very high energies {E^^ > 3 5 • 10^ GeV) in Fig. 3-5 
are explained by the insufficient simulation sampling for at- 
mospheric neutrino energy spectra (fig. 2) at corresponding 
energies and accuracies of multi-dimensional integrations. 
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